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The continuous progress in fabricating low-dimensional systems with large spin-orbit couplings 
has reached a point in which nowadays materials may display spin-orbit splitting energies ranging 
from a few to hundreds of meV. This situation calls for a better understanding of the interplay 
between the spin-orbit coupling and other interactions ubiquitously present in solids, in particular 
when the spin-orbit splitting is comparable in magnitude with characteristic energy scales such as 
the Fermi energy and the phonon frequency. 

In this article, the two-dimensional Frohlich electron-phonon problem is reformulated by introduc- 
ing the coupling to a spin-orbit Rashba potential, allowing for a description of the spin-orbit effects 
on the electron-phonon interaction. The ground state of the resulting Frohlich-Rashba polaron is 
studied in the weak and strong coupling limits of the electron-phonon interaction for arbitrary val- 
ues of the spin-orbit splitting. The weak coupling case is studied within the Rayleigh-Schrodinger 
perturbation theory, while the strong-coupling electron-phonon regime is investigated by means of 
variational polaron wave functions in the adiabatic limit. It is found that, for both weak and strong 
coupling polarons, the ground state energy is systematically lowered by the spin-orbit interaction, 
indicating that the polaronic character is strengthened by the Rashba coupling. It is also shown 
that, consistently with the lowering of the ground state, the polaron effective mass is enhanced 
compared to the zero spin-orbit limit. Finally, it is argued that the crossover between weakly and 
strongly coupled polarons can be shifted by the spin-orbit interaction. 

PACS numbers: 71.38.-k, 71.38. Fp, 71.70.Ej 



I. INTRODUCTION 

The Frohlich Hamiltonian describing a single electron 
coupled to longitudinal optical phonons is a paradig- 
matic model of the electron-phonon (el-ph) interaction^ 
and has represented in the past, in addition to its in- 
terest for the solid-state physics, an ideal problem for 
testing many mathematical methods in quantum field 
theoryi^ Because of the coupling with the phonon field, 
the resulting quasi-particle, the polaron, has an effec- 
tive mass larger, and a ground state energy lower, than 
the free electron. These quantities have been investi- 
gated for the three-dimensional (3D) case by means of 
perturbation theory for the weak-coupling limit p and of 
variational treatments for the intermediate;^ and strong- 
coupling casesi^ The path-integral variational calcula- 
tions of Feynman^ and subsequent refinements of this 
method^ have provided a solid description for all val- 
ues of the coupling, verified also by improved variational 
methods^ and by quantum Monte-Carlo studiesJ^— 

The interest aroused some time ago on semiconduc- 
tor heterojunctions, or other low-dimensional systems, 
prompted to modify the Frohlich model to accounting for 
two-dimensional (2D) and quasi-2D systemsii^ By apply- 
ing the same methods derived for the 3D case, the ground 
state properties for the strictly 2D case were evaluated for 
weak, strong and intermediate couplingS )^^d^'^^d^ and 
the obtained systematic lowering of the ground state en- 
ergy and the enhancing of the effective mass compared 
to the 3D case has pointed out the role of dimensionality 
in enhancing the polaronic character 



Concerning the el-ph problem in low dimensions, 
recent progresses in developing high-quality low- 
dimensional systems and in material engineering provide 
hints that, for a vast class of low-dimensional materi- 
als, the usual 2D Frohlich model, as considered in lit- 
erature, may be incomplete. This concern comes about 
by considering 2D systems exhibiting strong spin-orbit 
(SO) splitting of the electronic states due to the inver- 
sion asymmetry in the direction orthogonal to the con- 
ducting plane (Rashba SO mechanism). This situation is 
encountered in semiconductor quantum wells with asym- 
metric confining potentials^^ in the surface states of met- 
als and semi-metals,— i^i^i and in surface alloys such as 
Li/W(110),^ Pb/Ag(lll), .23,24 g^^^ Bi/Ag(lll),.^ with 
SO splitting energies ranging from a few meV in GaAs 
quantum wells to about 0.2 eV in Bi/Ag(lll)i^ In these 
systems, therefore, the SO energy may be of the same 
order or even much larger than the typical phonon fre- 
quency, rising the question of how such state of affair 
affects the el-ph interaction, in general, and the Frohlich 
coupling, in particular. 

As pointed out in several workS )^^i^'''^^i^^'^°i'^^ the 
Rashba interaction describing the SO coupling can have 
profound effects on the low energy properties of the itin- 
erant electrons. Namely, in the low density regime, the 
Rashba SO coupling induces a topological change of the 
Fermi surface of the free electrons, leading to an effective 
reduction of the dimensionality in the electronic density 
of states (DOS). In this situation, a 2D low density elec- 
tron gas would develop, in the presence of SO Rashba 
coupling, a phenomenology similar to one-dimensional 
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(ID) systems, triggered by the square-root divergence of 
the (effectively ID) DOS at low energies. 

Some interesting consequences of this scenario on 
the el-ph problem have already been discussed in 
Ref. |3C|. concerning the superconducting transition, and 
in Ref. [3l|] for the effective mass and the spectral proper- 
ties. The picture arising from these works, although be- 
ing limited to the momentum-independent Holstein el-ph 
interaction and to weak-to-moderate couplings, confirms 
that, for sufficiently low electron densities, the coupling 
to the phonons is amplified by the SO interaction through 
the ID-like divergence of the DOS. 

Notwithstanding the relevance of these results for the 
Holstein model, the use of a local el-ph interaction may 
however result inadequate in the extremely low elec- 
tron density regime, where the SO effects are more 
evident.— 1^ Indeed, the lack of effective screening in this 
case would rather suggest a long-range interaction as be- 
ing a more appropriate description of the cl-ph coupling. 
It is therefore natural to consider the 2D Frohlich po- 
laron, and its coupling to the SO interaction, as a model 
better describing the unscreened el-ph interaction in 2D 
Rashba systems in the low density limit. 

In this article, a single electron moving with a parabolic 
dispersion in the two-dimensional x-y plane is coupled 
simultaneously to the Rashba SO potential and to the 
phonon degrees of freedom through a Frohlich interac- 
tion term. The total system is then described by the 2D 
Frohlich- Rashba Hamiltonian H = H^i + Hph + Hei-ph, 
where {h — 1) 



with 



Hei^^ + 0(p) • a- 
Zm 



(1) 



is the Hamiltonian for an electron with mass m and mo- 
mentum operator p = — iV with components {px,Py,0), 
cr is the spin-vector operator with components given by 
the Pauli matrices, and r2(p) is the SO vector field which 
in the case of Rashba coupling reduces to: 



n(p) 




(2) 



where 7 is the SO coupling parameter. The phonon part 
of the Hamiltonian is given by 



Hph = '^o ^ aqflq, 



(3) 



where (aq) is the creation (annihilation) operator for 
a phonon with momentum q = {cix,qy) and optical fre- 
quency LjQ. The el-ph interaction Hamiltonian for the 
2D electron coupled to longitudinal optical (LO) phonons 
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Mo = it^o 
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where a = e^{e^ — Eq ) ■\/m/2a;o is the dimcnsionless el- 
ph coupling constant, with e being the electron charge, 
and £00 and eo the high frequency and static dielectric 
constants, respectively. 

It is worth clarifying here the significance of the 2D 
Frohlich interaction of Eq.([3]) with respect to the char- 
acteristics of specific materials. For quantum wells and 
2D heterostructures, where the electron wave function is 
assumed here to be confined in a sheet of zero thickness, 
Eq. (m describes the couphng of the electron to bulk LO 
phonons, while the coupling to interface phonon modes is 
neglected. The inclusion of such interface phonon contri- 
butions may be important in describing specific materi- 
als, but it is unnecessary for the present study, where the 
focus is on the SO effects on the unscreened (long-range) 
el-ph interaction, for which Eq.(|3]) is a paradigm for the 
2D case. Concerning the el-ph coupling of electronic sur- 
face states, Eq.Q coincides (apart from a redefinition 
of Mo) with the coupling to 2D surface phonons when 
the coupling to bulk phonons extending below the sur- 
face is negligible Such approximation is coherent with 
the ideal 2D assumption for the electron wave function, 
which is physically realized when the electronic surface 
states have negligible coupling to the bulk. A further mo- 
tivation of using the 2D Frohlich model ^ is that, in the 
absence of SO interaction, the ground state polaron en- 
ergy Ep and effective mass m* have already been studied 
by several authorS j^^i^^'^'^i^^'^^ and the exact results ob- 
tained for the weak {a <C 1) and strong (a ^ 1) coupling 
limits provide a useful reference for the effect of nonzero 
SO coupling. 

In the present work, the 2D Frohlich-Rashba Hamil- 
tonian is studied by considering the weak and strong 
coupling limits of the el-ph interaction, with arbitrary 
strength of the SO coupling 7. For a ^ 1 the polaron 
energy Ep and the effective mass m* are obtained from 
second order perturbation theory in SeclHl where nu- 
merical and exact analytical results are presented. It is 
shown that the effect of 7 ^ is qualitatively similar to 
that observed in the Holstein modelf^i^ namely, the SO 
coupling enhances the effective coupling to the phonons. 
In particular, Ep is lowered by 7 and, simultaneously, the 
effective mass m* is enhanced. In Sec lIIII thc strong cou- 
pling limit a ^ f is treated by the variational method, 
providing a rigorous upper bound of the ground state 
energy for arbitrary values of the SO interaction. As for 
the weak el-ph coupling case, it is found that Ep (m*) is 
lowered (enhanced) by the SO interaction, implying that 
the Rashba coupling always amplifies the polaronic char- 
acter, regardless of whether the el-ph interaction is weak 
or strong. 
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II. WEAK COUPLING 

111 the presence of SO interaction, the electron wave 
function is a spinor and its Green's function is conve- 
niently represented by a 2 x 2 matrix in the spin sub- 
space. For a = the free electron Green's function Go 
is readily obtained from H^i: 



Go(k,w) = [uj-- — n{k)-(T 



lY.[l + sn{k)-cT]Gt>{k,u), (6) 



s=± 



where k is a 2D electron wavenumber, f2(k) 
f2(k)/|f2(k)| and 



1 



UJ — kP' /2m — s^k 



(7) 



is the free electron propagator for the two (s — ±1) chiral 
states characterized by two distinct bands with shifted 
parabolic dispersions k"^ /2m ± ^k. The lowest band has 
its minimum value —Eq at k ^ kQ, where fco a-nd Eq are 
the Rashba momentum and energy defined respectively 
by: 



fco = rnj, Ef) = —7^ 



(8) 



For later convenience, it is useful to express the electron 
energy relative to Eq , so that the poles of Eq. ((7| appear 
at energies: 



£;±(fc) = -L(fc±fc„)2. 

Zm 



(9) 



The free electron ground state is then given by the elec- 
tron occupying the lower band at wavenumber k ~ ko 
with energy uj = 0. 

In the weak el-ph coupling limit (a ^ 1) the ground 
state properties are obtained by the electron self-energy 
evaluated in the second order perturbation theory. At 
zero temperature, the resulting single electron self-energy 
is therefore: 



S(k,tj) = \Mo 



(27r)2 |k-k 



-Go(k',c^-c^o)- (10) 



Because of the momentum dependence of the Frohlich 
interaction, and contrary to the Holstein el-ph case con- 
sidered in Ref . [Slj , the self-energy is not diagonal in the 
spin subspace. However, since the momentum depen- 
dence enter only through the modulus of the momentum 
transfer, equation (jlOp can be rewritten in a quite simple 
form. By using (0(k) • cr)^ = 1 and 

(fJ(k) • cr)(fj(k') • cr) = k • k' + (k X il)a^ay, (11) 

then the quantity Jl(k') - cr appearing in Eq. (fTU|) through 
Go(k', UJ — loq) can be replaced simply by (J^(k) • cr) k • k' 



because the second term of Ea. (|TT|) vanishes after the 
integration over k'. In this way, the resulting self-energy 
reduces to: 

S(k,u;) =Ed(fc,t^)l + So(fc,t^)f^(k)-cr, (12) 

where 1 is the unit matrix and and So are, respec- 
tively, the diagonal and off-diagonal contributions to the 
self-energy, both depending solely on the modulus of kj^ 
Their explicit expressions are: 



Sd(fc,t^) = 



So(fc,a;) = 



' f dk' ^ 1 
J (27r)2 |k-k' 



1 



|A/oP f dk' 



^ llr - W' 



\LU-LUO-E,{k')' 

(13) 

sk - k' 



2 J {2t:)^ ^ \k-k'\Lu -LUo-Es{k') 



(14) 



In the limit of zero SO coupling, since Es{k) k^/2m, 
So(fc,a)) vanishes because of the summation over s = 
±1 in Eq. (fTi|) . Notice also that, independently of 7, 
So(fc,a;) = when the factor l/|k - k'| in Eq.lUll is re- 
placed by a constant, as in the momentum- independent 
Holstein el-ph coupling model. 

By using Eg. p^ the Dyson equation for the interact- 
ing propagator G reduces to 

G"i(k,cj) = Go\k,uj) --Eik, uj) 

^ UJ - Sd(fc,Cj) - Eq 

Zm 

-[jk + j:o{k,uj)]n{k)-(T, (15) 

and the poles uj± of G are then given by: 

uj±^ E±{k) + j:4k,uj±)±^oik,uj±). (16) 

Now, the Rayleigh-Schrodinger perturbation theory per- 
mits to evaluate the lower energy pole uj- at the lowest 
order in the el-ph coupling a. This is accomplished by 
replacing uj- by the unperturbed energy E-{k) in the 
energy variables of and So- In this way, the lower 
pole reduces to w_ = E^{k) + S_(A:) -I- 0(0^), where 

S_ (fc) = Sd(fc, S_ (fc)) - So(fc, E^ (fc)). (17) 

Finally, by expanding S_(fc) up to the second order in 
k~ ko, the polaron dispersion in the vicinity of ko can be 
written as: 



uj.^Ep + ^ik-k*f, 



(18) 



where the polaron ground-state energy Ep, the effective 
mass m*, and the effective Rashba momentum /cq are 
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jiven respectively by: 



the effective mass can be expressed as 



m 
m 



"-0 
fco 



S_(fco)-^S'_(fco)' 
E_(fco)+O(a2), 

[l + mSl(fcor' 

1 -mS':(fco) + 0(a2), 

l-^I]'_(fco) 

l-^S'_(fco) + 0(a'). 

fco 



(19) 
(20) 

(21) 



Let us first consider Ep and m* . In the zero SO limit, 
Eqs. HH) and (^0]) at fco = lead respectively to iJp = 
7raa;o/2 and m* /m = 1 + 7ra/8, which correspond to 
the results already reported in Refs. ^13, 14,15]. For finite 
values of the SO coupling the ground state energy and 
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FIG. 1: (a): ground state energy factor fEpie) as a function 
of the SO parameter eo = Eq/ujo. The solid line is the nu- 
merical calculation, while the dashed line is the weak SO limit 
Ea. (|25|l . Inset: fspie) plotted for a wider range of eo- (b): 
the effective mass factor /m*(eo) from numerical integration 
(solid line) and from Ea. (|26|l (dashed line). Inset: f,n*{eo) 
plotted for a wider range of £o- 



Ep = --aujofEpieo), 

1 + -a/„i.(eo), 

m 8 



(22) 
(23) 



where the factors /Ep(eo) and /m«-(eo) contain all the 
effects of the SO interaction and depend solely on the 
dimensionless SO parameter 



Eo mj^ 
eo = — -z — 

Wo ^Wo 



(24) 



In the weak SO limit, the self-energy terms ([T^ and (fTil) 
can be expanded in powers of the SO interaction, allowing 
for an analytical evaluation of the integrals. In this way, 
up to the linear order in Eq, /_Ep(£o) ^^nd fm'i^o) are 
found to be: 



/£p(£o) = l + f +O(£o), 



(25) 
(26) 



indicating that the polaronic character is strengthened 
by the SO interaction since, through Eqs. ([22|) and (|23p. 
the polaron energy Ep is lowered and, simultaneously, 
the effective mass m* is enhanced when Sq > 0. This fea- 
ture is not limited to the small eo limit, but holds true 
for arbitrary strengths of the SO coupling. This is shown 
in Fig. [1] where /ep(£o) and /m*(eo), obtained from a 
numerical integration of Eas. p3|) and (|14p . are plotted 
as a function of Eq by solid lines and compared with Eqs. 
(PS)) and (pS)) (dashed lines). The same quantities calcu- 
lated for a wider range of £o are plotted in the insets of 
FigO and confirm that the ground state energy Ep and 
the effective mass m* are continuous functions of eo and 
are, respectively, further lowered and enhanced by the 
SO couphng. In the strong SO hmit £o ^ 1, it is found 
that fEri^o) grows as In(eo) while /„• (eq) grows linearly. 
It is interesting to note that the Holstein-Rashba model 
studied in Ref. [sij predicts results qualitatively similar 
to the Frohlich model, indicating that the SO interaction 
strengthen the polaronic character independently of the 
specific form of the el-ph interaction!^ 

In addition to Ep and m*, the interplay between the 
el-ph coupling and the SO interaction modifies also the 
Rashba momentum fco through Eq. (j2ip . In the weak SO 
limit, the effective quantity fcQ is found to be 



fco 



32 



(27) 



indicating a reduction of the bare Rashba momentum fco, 
confirmed also by the numerical calculation of Eq. (PT|) 
reported in Fig. [2] by the solid line. As shown in the 
inset, for fixed el-ph coupling a, fcg however does not 
deviate much from its bare limit fco, even for large values 
of the SO parameter eo. 

Let us compare now the present results with those 
appeared recently in literature. In Ref. [35| the ground 
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FIG. 2: Effective Rashba momentum kg as a function of the 
SO parameter eo = Eq/ujq. The numerical integration of 
Ea. (|21|l (soUd line) is compared with the weak SO result (|27p 
(dashed line). Inset: the same quantity plotted for a wider 
range of eo- 



state energy of a polaron near a polar-polar semicondue- 
tor interface with Rashba SO coupling has been evalu- 
ated with the Lee-Low-Pines method*i As a function of 
the SO splitting, the polaron ground state is found to 
be lowered, in qualitative agreement therefore with the 
present results. A more quantitative comparison is how- 
ever precluded by the different model of Ref. [s^ , where 
contributions from interface phonon modes and confin- 
ing potentials are considered as well. In another work^ 
the Rayleigh-Schrodinger perturbation theory has been 
applied to the polaron ground state of the 2D Frohlich- 
Rashba model, permitting therefore a direct comparison 
with the analysis presented here. Despite that the au- 
thors of Ref. [301 find that the polaron ground state is 
lowered by Eq, their values of Ep differ from those plot- 
ted in Fig. [TJa). In Ref.[3^, in fact, the ground state 
energy factor is found to be fsp i^o) = — £o, 

which implies a small Sq expansion different from Eq. (|25p 
and, more importantly, a divergence of Ep at eg = 1- In 
Fig. lUa), instead, nothing of special happens at eo = 1- 
This discrepancy is easily traced back in the fact that in 
Ref. [1^ the expansion of S_(fc), Eq. pT]) . is made around 
fc = 0, instead of fc = fco as done here, which does not cor- 
respond to a pcrturbativc calculation of the ground-state 
energy. 

The results presented in this section have been derived 
by assuming a weak coupling to the phonons. However, 
as it is clear from the plots in Fig. [1] the enhancement of 
the polaronic character driven by Sq for fixed a unavoid- 
ably renders the perturbative approach invalid for suffi- 
ciently large Eq values. For example, from Eq. (P5)l . the 
validity of the weak coupling results for m* /m are sub- 
jected to the condition af„i'{eo) ^ 1, otherwise higher 
order el-ph contributions should be considered for a con- 
sistent description of the SO effects. The question re- 
mains therefore whether the SO enhancement of the po- 



laronic character survives also for large a values, or it is 
instead limited to the weak coupling limit. In the next 
section, this problem is studied for the limiting case of 
strong cl-ph interaction a 3> 1, providing therefore, to- 
gether with the weak coupling results, a global under- 
standing of the SO effects on the Frohlich polaron. 



III. STRONG COUPLING 

It is well known that a perturbative scheme such that 
employed in the previous section fails to describe the 
Frohlich polaron ground state when the el-ph coupling 
is very large. This is due to the fact that for a ^ 1 
the lattice polarization, and resulting "self-trapping" ef- 
fect experienced by the the electrouf^ renders the plane 
wave representation of the unperturbed electron inappro- 
priate for obtaining the polaron ground state. Instead, 
as originally proposed in Ref. @ and rigorously proved in 
Refs. [38ll39| . the asymptotic description of the polaron 
wave function in the strong coupling limit a 3> 1 is 
that of a product between purely electronic, ■(/'(r), and 
purely phononic, |^), wave functions. Within such adi- 
abatic limit, the ground state energy and the effective 
mass of a 2D Frohlich polaron have been calculated in 
Refs. [11113 by using the variational method with differ- 
ent ansatz wave functions. From Ref. [l^ . one realizes 
that exponential, gaussian and Pekar-type wave functions 
provide increasingly better estimates of Ep with accura- 
cies respectively of 14%, 0.3%, and 0.03% with respect to 
the exact ground state energy Ep/luq — — 0.40474q;^, ob- 
tained by a numerical solution of the integro-differential 
equation for the electron wave functioUfi^ In the follow- 
ing, the variational method is used to evaluate the SO 
effects on the polaron ground state. 



A. trial wave functions 

For the nonzero SO case, due to the presence of the 
Pauli matrices in Eq.([T]), suitable ansatz wave functions 
must take into account the electron spin degrees of free- 
dom. Hence, in full generality, the strong-coupling po- 
laron wave function may be represented as: |'4',0 = 
^(r)l^), where ^'(r) is a two-components spinor for the 
electron. The corresponding expectation value of the to- 
tal Hamiltonian H is: 



(*,^ii?i*,c) =(*|i/,;i*) + (eiifp,ie) 



+ ^T. ^(*^op(q)(Cl«qlO + h-c), 

(28) 



where 



p{q) = (^le'i-n*) = / fire''i"'|*(r)| 



(29) 
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The form of Eq. ([28)1 permits to integrate out the phonon 
wave function in the usual way. Hence, by introducing 

the phonon coherent state |^) — A/'e^<i^'''''i|0), where M 
is a normalization factor and a variational parame- 
ter, minimization of (|28p with respect to leads to the 
functional 

Em = p^-\pm\ (30) 

where the continuum limit A~^J2ci ^ /'^'^/C^"')^ 
been performed. By choosing an appropriate functional 
form for ^'(r), and by minimizing E[^] with respect 
to the variational parameters defining 'S'(r), an upper 
bound for the ground state energy is then ii'[\I'o], where 
\I'o(i") is such that i?[*o] = min(i?['4']). As done in 
the previous section, the polaron energy is then obtained 
from 

Ep=E[<^o]+Eo, (31) 

where Eq is the free-electron SO energy defined in Eq.®. 

Of course, the functional form of ^(r) is decisive for 
obtaining accurate estimates of the ground state energy, 
and a suitable choice must be guided by looking at the 
properties of the true ground state spinor ^0(1"). These 
can be deduced by a formal minimization of the func- 
tional E[^] with respect to ^. By introducing the La- 
grange multiplier e to ensure that the wave function is 
normalized to unity, minimization of ([50]) leads to: 

i?ei*(r) + Vir)^{v) = e*(r), (32) 

where, by using the definition of p(q) given in in Ea. (|29p : 

_ _2Wol f dq P(q)* .... 

= .\M£f,,\p^. (33) 
TTUJo J |r-r'| 

I 



where U{r) = V{r) — V/2 and the polaron energy is ob- 
tained from Ep = Eq + Eq. By introducing the dimen- 
sionless variable p = r/ip, where £p = l/a(mwo)"'^/^ is a 
measure of the polaron spatial extension in the zero SO 
limit, and by noticing that Eq does not depend on the 
sign of 7, it is straightforward to realize from Eqs. ([36)) 



From the above expression of ^(r), the functional pO)) 
can be rewritten as E[^] = (*|ifez|*) -I- V/2, where 
V = (*|y(r)|*). Now, if *G is the exact ground state 
wave function, with ground state energy Eq = £'['S'g], 
then, from (l32l)_and Eg = + V/2, it is found 

that e = Eg + V/2, so that Eq.§^ reduces to: 



i7e/*G(r) + [V{r) - y/2]*G(r) = SG*G(r). (34) 



As noted in Ref.flil (see also Refs.jHiH), the ground- 
state wave function of a 2D electron subjected to a SO 
Rashba interaction and to a 2D central potential (i.e. a 
potential depending only upon r ~ |r|) is of the form 



where is the azimuthal angle of r. Now, if Eq. pS]) 
is used in Eq.(|33|). the resulting self-consistent potential 
depends only upon r, V[v) V{r), so that Eq. ((35|) is 
consistently also the correct form for the polaron ground- 
state wave function. Hence, passing to polar coordi- 
nates, Eq. (p4|) can be rewritten as a system of integro- 
differential equations for the spinor components ipi and 

1p2- 



(36) 
(37) 



and (j37p that the polaron ground state energy scales as 

Ep=T(^^)a^iUo, (38) 

where Eq = Eq/ujq is the dimensionless SO energy intro- 
duced in Eq. ((24|) and is a generic function. It is found 
therefore from Eq. ([38|) that the dependence of Ep on the 
SO interaction is through the effective parameter Eq/o? , 



1 



1 d 



2m \ dr'^ ^ r dr ^ 



V'i('') - 7 ( T" + - ) '^2{r) = EGipi{r), 



1 



1 d 



2m V c?r^ ^ r dr 1 ^ ^ 



dr 



V'2(f) + l-r^^i{r) 
dr 



EG'ip2ir), 
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which is treated in the following as an independent vari- 
able. Although Eq/o^ is then formally allowed to vary 
from to oo, it is nevertheless important to estimate the 
range over which e^/a^ is expected to vary for reason- 
able values of the microscopic parameters -Bq, ojq, and a. 
To this end, it must be reminded that the strong cou- 
pling limit of a 2D Frohlich polaron (in the absence of 
SO interaction) is appropriate only for a ^ 5?^^ and that 
the typical phonon energy scale is of the order of few to 
tens meV, say « 5 — 10 meV. The largest value of the 
Rashba energy f^o reported so far is of about 0.2 eV?^ 
so that £o/a^ ^ 1 — 2 is a rather conservative estimate 
compatible with material parameters and with the strong 
coupling polaron hypothesis. 

Let us now evaluate the behavior of ipi (r) and V'2 (f ) for 
r <^ £p and r ^ £p. By requiring a regular solution at 
the origin, it turns out by inspection of Eqs. ((36|) and (p7|) 
that the spinor components of (|35p behave as ipi{r) = 
const, and tpiii^) oc r as r ^ 0, while the behavior for 
r ^ is obtained from the large r limit of Eqs. ([5B|) and 



1 d^Mr) yMr) 
dr 
dipiir) 



2m dr^ 

1 d'^i!2{r) 
2m dr^ 



1 



+ 7 



dr 



(39) 
(40) 



where the quantity W = Eq + V /2 is negative for bound 
states. Solutions of Eas. p9|) and (|40p which are finite 
for r — > 00 are linear combination of cxp(— A+r) and 
exp(— A_r) with 



A+ = 



-2m{Ep + V 12) ± i/fco 



(41) 



implying an exponential decay of the polaron wave func- 
tion, accompanied by periodic oscillations of wavelength 
27r/fco. 

The informations gathered on the limiting behaviors 
of the ground state wave function are sufficient for guess- 
ing some appropriate trial wave functions to be used in 
Eg. ([50)1 . By assuming that for zero SO coupling the elec- 
tron is in a spin-up state, then a simple ansatz compatible 
with the limits discussed above is 



*(r) = /(r) 



cos(6?-) 
sin(6r) e*"^ 



(42) 



where 6 is a variational SO parameter vanishing for 7 = 
and /(r) is an exponentially decaying function for r ^ 00 
and such that /(O) 7^ 0. The advantage of Eq.fgll) is that 
one can use exponential or Pekar-typc functions for /(r). 
automatically recovering therefore the known results for 
the zero SO case>i^ It should be noted, however, that 
in the U(r) — *■ limit Eq. (|l^ does not reproduce cor- 
rectly the behavior of the exact ground state wave func- 
tion, which is instead given by Ea. (|35p with "01 (r) and 
'^i{r) proportional to the Bessel functions Jo{kor) and 
Ji(fco'')j respectively4Si^ Hence, Ea. (|^ is not expected 
to provide a reliable ground state energy in the strong SO 



regime, for which U{r) can be treated as a perturbation. 
To remedy to this deficiency, the following alternative 
form of the polaron ansatz is proposed: 



*(r) - /(r) 



Mbr) 
Ji{br)eJ'f 



(43) 



where, as before, 6 is a variational SO parameter. As it 
will be shown below, the lowest value of Ep is given either 
by Eq. (|^^ or by Eq. ([151) ' depending on the specific form 
considered for /(r) and on the value of the SO coupling. 



B. ground state energy 

To evaluate the polaron ground state energy, three dif- 
ferent trial wave functions for f(r) are considered: ex- 
ponential, Gaussian and Pckar-type. As shown below, 
the Gaussian ansatz will provide results comparable to 
those coming from the exponential and Pekar functions, 
despite its faster decay for r 00 compared to Eq. (|¥T|) . 
These three trial wave functions will be used in combi- 
nation with the sinuisodal and the Bessel-type spinors 
of Eqs. (|42l) and (|43)) . respectively, giving a total of six 
different ansatzes for the Frohlich-Rashba polaron wave 
function. 

Exponential ansatz. Let us start by evaluating the 
functional Eq. pO)) . by using the exponential ansatz 

f{r) = ^exp(— ar), where a is a variational parameter 
and ^ is a normalization factor, in combination with the 
sinuisodal trial wave function (|42p . By introducing the 
dimcnsionless quantities d = alp,b ~ Mp, and 7 = k^lp, 
for nonzero SO interaction the functional (|30p evaluated 
with the exponential-sinuisodal ansatz reduces to 




(44) 



For weak SO couplings, Eq. (|44p has its minimum at 6 = 
7 = \J2ei^l a and d ~ 3\/27r/16, so that the resulting 
polaron energy Ep = £'['4'o] + E^ becomes 



Ep 



£0 



o 



(45) 



In the £0 = limit, Eq. (|i5)) reduces to Epjop-uj^ = 
— (37r/16)^ — —0.3469, recovering therefore the result of 
Ref. [iJI, while for £0 > the polaron energy is lowered 
by the SO interaction, in qualitative analogy with the 
weak electron-phonon behavior discussed in Sec [ill The 
lowering of Ep is confirmed by a numerical minimization 
of Eq. fi^ whose results are plotted in Figini^a) (open cir- 
cles). For £o/a^ = 1, the polaron energy has dropped to 
EpjcP'iijQ ~ —0.65, that is about two times lower than 
the zero SO case. However, upon increasing Sq/o^, Ep 
displays a minimum at Eq/o^ ~ 3.98 [inset of FiglJl^a)] 
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FIG. 3: Polaron ground state energy as a function of eo/o? for different trial wave functions for /(r). (a): exponential; (b): 
Gaussian; (c): Pekar. The sinuisodal and the Bessel type of anstazes are given respectively by Eq.(j42} and Ea. (|43|l . Inset: the 
polaron energy for a wider range of £0/0? values. 



and for larger values of the SO interaction the polaron 
energy increases. Eventually, for £q/o? > 14 the calcu- 
lated ground state energy becomes larger than the zero 
SO value Ep/a^ujo = — (37r/16)^. Such upturn of Ep 
for large stems from the inadequacy of the sinuiso- 
dal components of (j42|) in treating the oscillatory be- 
havior in the strong SO regime which, as pointed out 
above, should instead be given by Bessel-type functions. 
Indeed when the exponential ansatz for /(r) is used in 
Eq. (j43|) . rather than in Eg. (|42|) . not only the resulting 
Ep is lower than the previous case, but also the upturn 
of Ep disappears, leading to a monotonous lowering of 
the polaron energy as s^jo? increases [filled circles in 



FigUa)]. As £o/c 



00, however, the polaron energy 



does not decrease indefinitely but rather approaches a 
limiting value. Although an accurate numerical evalua- 
tion of Ep for EqIo? > 100 has turned out to be difficult, 
the asymptotic value of Ep can nevertheless be obtained 
analytically from the strong SO limit of the exponential- 
Bessel expression for 



(46) 



whose minimum is at 5 = 7 and a = n/ a/2, leading to 
Ep 



lim 



~ -2.467. 



Gaussian ansatz. The results obtained by using a 
Gaussian wave function of the form /(r) — ylexp(— a^r^) 
are plotted in Fig. Elb). Compared to the exponential 



wave function, the Gaussian ansatz gives an overall low- 
ering of the polaron energy for both sinuisodal and Bessel 
forms of the spinors. In the Eo/a^ <C 1 limit, and inde- 
pendently of which particular spinor is used, the ground 
state polaron energy is found to be: 



Ep 



£o_ 

r.2 



O 



(48) 



confirming in this regime the linear dependence on the SO 
coupling of Eq. (|l5)) . For larger values of the SO coupling, 
and contrary to the case shown in Fig. OJa), the sinuiso- 
dal and Bessel-type spinors give basically the same values 
of Ep for all SO couplings up to E^jo? ~ 1. Beyond this 
value, as for the case with the exponential wave function, 
the polaron energy obtained from the sinuisodal ansatz 
becomes larger than that obtained from the Bessel spinor 
and, as shown in the inset of Fig. EJb), rapidly increases 
while the Gaussian-Bessel anstaz gives a monotonous 
lowering of Ep. For E^jc? ^ 1, the Gaussian-Bessel 
energy functional has the same form of Eq. (|46)) with the 
latter term substituted by —2.279a, which implies 



lim 



Ep 



eo/ofl^oa a UJQ 



-2.579. 



(49) 



Pekar-type ansatz. Let us now evaluate Ep by using 
in Eqs.lHl) and ^ the Pekar-type ansatz /(r) = A{1 + 
air + a2r^) exp(— ar). For zero SO coupling, this ansatz 
gives Ep/a^uio — — 0.4046^^ which is a lower energy 
than those obtained from the exponential and Gaussian 
trial wave functions and only 0.03% higher than the exact 
result -0.40474 of Ref.fTl]. As shown in Fig. E^c), the 
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Pckar-type ansatz gives slightly better estimates of Ep 
also for nonzero SO couplings, with an overall behavior 
similar to the previous cases. Namely, in the weak SO 
regime one finds 



Ep 



-0.4046 



£o 



o 



(50) 



and, as before, for stronger SO couplings the energy ob- 
tained from the sinuisodal spinor increases indefinitely 
with eo/a^. However, contrary to the exponential and 
Gaussian ansatzes, the Pekar-type wave function may 
give a lower polaron energy when used in combination 
with the sinuisodal spinor. This holds true as long as 
Eo/a^ < 2.72, while for stronger SO couplings it is the 
Bessel-type spinor which gives the lower Ep [inset of Fig. 
[31(c)]. A numerical minimization of the asymptotic limit 
of the Pekar-Bessel functional for Eq jo? ^ 1 gives 



lim 



Ef 



-2.91, 



(51) 



which is lower than the asymptotic values of Eqs. (|T7|) and 
dMl). 

The results plotted in Fig. [3] clearly demonstrate that, 
since the variational method provides an upper bound for 
true ground state polaron energy, the lowering of Ep in- 
duced by the SO coupling is a robust feature of the strong 
coupling Frohlich-Rashba polaron. Among the different 
ansatzes studied, the lower polaron energy is obtained by 
using a Pekar-type wave function for /(r) in combination 
with the sinuisodal spinor for weak to moderate values of 
Eq/o^ or with the Bessel-type spinor for stronger SO cou- 
plings. Given that, as discussed above, reasonable values 
of e^ja^ for strongly-coupled polarons fall in the range 
< EqIo? < 1 — 2, the Pckar-sinuisodal wave function 
provides therefore the best description of the Frohlich- 
Rashba polaron in this regime. 



C. effective mass 

As demonstrated in Scc|TTl the effective mass m* of 
a weakly-coupled polaron is enhanced by the SO in- 
teraction and, given the results above, the same phe- 
nomenon is reasonably expected to occur also for the 
strong-coupling case. To quantify the polaron mass en- 
hancement within the localized wave function formal- 
ism, it is useful to follow the approach of Refs. 4^l43ll4^ . 
briefiy described below, where a moving wave packet is 
constructed from the localized wave function. The quan- 
tity to minimize is 



j„[*',r] = (*',rii?-'^-p|*',e') 



(52) 



where is a Lagrange multiplier, which will turn out to 
be the mean polaron velocity, and P = p + q aJjCq is 
the total momentum operator. The wave function |'4'', 
is given by the product ^^'(r)!^') where 

*'(r) = e*P«-""*(r) (53) 



is the electron wave packet with po being a variational 
momentum, 'S'(r) is the ansatz localized wave function, 

and 1^') = 7Ve^<i^^''^|0). Minimization of dS!]) with re- 
spect to gives now the functional 



dq |p(q)? 



1 



(27r)2 q cjo-q- (^Vl*')' 

(54) 



where p is the electron momentum operator and /o(q)' = 
(\I>'|e*i''|\I/'). By using Eq.flHSl), it is easily shown that 
Jt,['3/'] reduces to 



J.[*']=(*|ifeH*> + ||-PO-t^ 



I Mo I 



rfq |p(q)l' 



(27r)2 q tjo-q-1^' 



(55) 



where p(q) = i^\e^'^'^\^) . Equation is minimized 
with respect to po by setting po = mv and, by expanding 
the last term of Eg. ([55]) up to the second order in tj, the 
corresponding minimum Jtj[^] becomes!^ 



2|Mn 



<^q (q ■ u)' 

(2^)2 
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|p(q)l' 



(56) 

where E\^\ is given in Eg. pop . From the above expres- 
sion, it is clear that Jt,[>&] differs from Jq\^\ at least to 
order . Hence, if and 'S'o ai'c the wave functions 
which minimize and Jo[\I'], respectively, then the 

difference ^f^; — 'S'o is also of order . As a consequence, 
the minimum of ([5^)1 . differs from Ju['4'o] only to 

order — ^o)^ = 0{v^^ so that, by neglecting terms of 
higher order than u^, minimization of (j56p is achieved by 
the best wave function which minimizes i?['4']. Therefore, 
by using i?[*o] ^ Ep ~ Eq and evaluating (\I>o|P|*o), 
from Eas. (|52)) and (|56|) it turns out that 



Ep{v) = Ep- 



1- 



2|Mo|2 f dq (q-u)2 



(2^)2 



|po(q) 



(57) 

permitting us to identify the quantity within square 
brackets as the mass enhancement factor m*/m. By in- 
tegrating over the direction of q and by using ([5]), m* /m 
becomes in the strong-coupling limit 



/27ra 



m (mijJo)'^/'^ 



(58) 



which; by replacing the momentum variable by the di- 
mcnsionless quantity q ~ q£p, gives a mass enhancement 
proportional to in the zero SO case. By using the ex- 
ponential, Gaussian, and Pekar-type ansatzes in Eq. (|58p . 
the resulting mass enhancement factor becomes m* /m = 
(3/16)Ma4 ^ o.6421a'^, m* /m = (j^j^fa' ~ 0.6170^, 
and m* jm ~ 0.73a^, respectively.— 

The results for nonzero SO coupling are plotted in 
FigHl for the sinuisodal (open circles) and Bessel (filled 
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FIG. 4: Polaron mass enhancement m* jra in units of oi' as a function of eq/o? for different ansatz wave functions, (a): 
exponential; (b): Gaussian; (c): Pekar. Inset; ra* jma'' is plotted for a wider range of SO values. 



circles) spinors evaluated with exponential (a), Gaus- 
sian (b), and Pekar-type (c) wave functions. For all 
cases, m*/m increases with So/a^ without much quan- 
titative differences between the various ansatzes as long 
as Eq/o^ < 2. As shown in the insets of Fig. IH for 
larger values of the SO coupling the use of the sinuisodal 
spinor largely overestimates the increase of the effective 
mass compared to the Bessel-type spinor results. How- 
ever, despite of the weaker enhancement of m * /m, the 
Bessel-type spinors give nevertheless an infinite effective 
mass at St^jo? = oo. Indeed, independently of the par- 
ticular form of f{r), for e^jc? — > oo the expectation 
value ('5'o|e*'i""|'i'o) appearing in Eg. ([55)) goes hkc ajq 
for g ^ oo, rendering the integral over q of Eg. ((58)) di- 
vergent. 



IV. DISCUSSION AND CONCLUSIONS 

The results presented in the previous sections consis- 
tently show that, for both the weak and strong coupling 
limits of the el-ph interaction, the ground state energy 
Ep of the Frohlich-Rashba polaron is lowered by the 
SO interaction and the mass is enhanced, leading to 
the conclusion that the Rashba coupling amplifies the 
polaronic character. This scenario suggests also that a 
weak-coupling polaron at eg = may be turned into a 
strong-coupling one for > or, more generally, that 
the crossover between weakly and strongly coupled po- 
larons may be shifted by the SO interaction. This pos- 
sibility can be tested by looking at the curves plotted in 
the main panel of Fig. [51 where the weak and strong 
coupling results for Ep/loq are reported as a function of 



the el-ph coupling a for different Eq values. For Eq = Q, 
the polaron energy follows Ep/ujq ~ — 7ra/2 for small a 
and Ep/ivo ~ — 0.4046q;^ for large a. These two limiting 
behaviors arc plotted in Fig. [S] by the uppermost curves 
and compared with a numerical solutions of the Fcyn- 
man variational path integral for the 2D polaron (filled 
circles). The largest deviation of the path integral solu- 
tions from the weak and strong coupling approximations 
falls in the range of intermediate values of a and signals 
a region of crossover between the weakly and strongly 
coupled polaron. A rough estimate of the crossover po- 
sition is given by a "critical" coupling, say a*, obtained 
by equating the weak and strong coupling results. For 
eo = therefore one has 7ra/2 = 0.4046q!^, which gives 
a* ~ 3.9. Now, as shown in Fig. [Slfor Eq = 5 and Eq = 20, 
the increase of the SO interaction systematically reduces, 
for fixed a, the polaron ground state energy and, at the 
same time, shifts the intersection point between the weak 
and strong coupling curves towards smaller values of the 
el-ph interaction. The "critical" value a* of the crossover 
is therefore reduced by the SO interaction. For Eq — 5 
and £o = 20 it is found that a* ~ 3.6 and a* ~ 2.7, 
respectively. The systematic reduction of the crossover 
coupling by the SO interaction is made evident in the 
inset of Fig. [5l where a* is plotted as a function of Eq- 
From Fig. [5] it is also expected that, beside the reduc- 
tion of a* , the crossover region is likely to be narrowed by 
Eq. Indeed, the intersection between the weak and strong 
coupling solutions for eo = 20 is apparently smoother 
than the case for £o = 0, suggesting that the true ground 
state energy would deviate less, and in a narrower region 
around a*, from the weak and strong coupling solutions. 

The scenario illustrated above, and in particular the 
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FIG. 5: (Color online). Ground state polaron energy Ep as 
a function of the el-ph coupling a for different values of the 
dimensionless SO parameter £o = Eq/uiq. The straight lines 
at small a refer to the weak coupling results, while the curves 
at large a are the solution of the strong-coupling theory. The 
filled circles are the solution of the Feynman path integral 
ansatz (see text). The point of intersection between the weak 
and strong coupling curves is a measure of the crossover el-ph 
coupling Q*. Inset: a* is plotted as a function of eo. 

SO effect on the crossover coupling, may be verified 
by quantum Monte-Carlo calculations of the Frohlich- 
Rashba action or, more simply, by generalizing the Feyn- 
man ansatz for the retarded interaction to eo > The 
results presented here on the limiting cases a <C 1 and 
a 3> 1 may then serve as a reference for such more gen- 



eral calculations schemes for arbitrary values of the el-ph 
coupling and of the SO interaction. 

Let us discuss, before concluding, possible generaliza- 
tions of the Frohlich-Rashba model employed here and 
the consequences on the polaronic character. Let us re- 
mind that in Ref . [3l| it has been demonstrated that also 
for a momentum independent el-ph interaction model, 
the Rashba SO term leads to an effective enhancement 
of the el-ph coupling. The SO induced lowering of the po- 
laron ground state is therefore robust against the specific 
form of the el-ph interaction, so that a similar behavior 
is expected to occur also when considering the contribu- 
tions from interface or surface phonon modes. However, 
a different form of the SO interaction term may lead to 
a much weaker effect. Consider for example the situa- 
tion in which, in addition to the Rashba SO coupling, 
the system lacks also of bulk inversion symmetry, as in 
III-V semiconductor hcterostructures, leading to an ex- 
tra SO term of the Dresselhaus typej ^^i^^ When both SO 
contributions are present, the square root divergence of 
the DOS at the bottom of the band of the free electron 
disappears, and it is replaced by a weaker logarithmic 
divergence at higher energies. In this situation therefore, 
at least for weak el-ph couplings, the SO interaction is 
expected to have a weaker effect on the polaron ground 
state, which tends to vanish as the Dresselhaus term be- 
comes comparable to the Rashba one. 

Let us conclude by noticing that, recently, the possi- 
bility of varying the coupling of 2D Frohlich polarons in 
a controlled way has been experimentally demonstrated 
by acting on the dielectric polarizability of organic field- 
effect transistors 4i The results presented here suggest 
that tunable 2D Frohlich polarons may be achieved also 
by acting on the SO coupling, which can be tuned by 
applied gate voltages in quasi-2D structured materials. 
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